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SUMMARY 

Bromodichloromethane  (BDCM)  was  tested  in  the  medaka 
fish  ( Oryzias  latipes)  model  in  a  short-term  and  a  long-term 
exposure.  The  short-term  exposure  lasted  <9  hours  with  the 
fish  tissue  and  blood  assayed  for  BDCM  at  the  conclusion  of 
the  exposure  and  this  data  was  used  to  build  a 
physiologically  based  pharmacokinetic  model  (PBPK)  for  BDCM 
uptake  in  medaka.  The  long-term  exposure  lasted  six  months 
and  included  a  48-hour  initiation/promotion  stage  with 
diethylnitrosamine  (DEN) .  At  the  conclusion  of  the  long¬ 
term  exposure,  the  fish  tissues  were  assessed 
histopathologically  for  determination  of  chemical  dosing 
effects . 

The  short-term  exposure  (USACEHR  Test  108-001)  used  the 
following  BDCM  test  concentrations:  0,  1.5,  and  15  mg/L. 
Twenty  same  gender  fish  were  dosed  in  each  1.5  L  of  test 
solution.  There  were  three  replicates  of  each  concentration 
and  gender.  Six-month  old  fish  were  used  for  this 
experiment.  Fish  length  and  weights  were  determined  at 
animal  necropsy,  as  well  as  weights  for  blood,  liver, 
spleen,  kidney,  gill,  and  muscle.  Samples  were  sealed  in 
glass  headspace  vials  and  then  analyzed  for  BDCM  content. 
Summarized  BDCM  tissue  data  was  provided  to  statistical 
modelers  and  a  PBPK  model  was  developed.  This  model 
describes  the  absorption,  disposition,  metabolism,  and 
excretion  of  BDCM  in  medaka.  Kinetic  predictions  from  the 
PBPK  model  were  linked  with  selected  pharmacodynamic 
responses . 

The  long-term  exposure  (USACEHR  Test  108-002)  tested 
the  following  BDCM  concentrations:  0,  1.5,  and  15  mg/L. 
Juvenile  fish  were  dosed  with  either  0  or  10  mg/L  of  DEN  for 
48  hours,  given  24  h  to  recover  in  processed  well  water, 
then  moved  to  the  BDCM  dosing  tanks.  Fish  assignment  was 
determined  by  DEN  initiation  and  by  BDCM  level,  with  three 
replicates  of  each  combination.  The  initial  number  of  fish 
in  each  flow-through  aquarium  was  30.  After  72  hours  of 
BDCM  dosing,  5  fish  were  removed  from  each  replicate  for 
labeling  with  5' -bromodeoxyuridine  (BrdU)  to  assess  cell 
proliferation.  The  selected  fish  underwent  a  72-hour 
exposure  to  75  mg/L  BrdU,  followed  by  euthanasia, 
sectioning,  staining,  and  counting  of  the  proliferating 
liver  cells.  Results  were  analyzed  statistically.  No  BDCM 
treatment  effects  were  found  for  cell  proliferation.  The 
remaining  25  fish  per  replicate  were  continuously  dosed  with 
BDCM  for  six  months.  Each  entire  fish  was  preserved  and 
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assessed  histopathologically  for  tissue  changes  between 
control  and  BDCM-dosed  groups.  Findings  were  analyzed 
statistically.  Briefly,  statistically  significant  findings 
were  found  in  the  15  mg/L  BDCM  groups,  regardless  of  DEN 
concentration.  These  significant  findings  included  changes 
in  the  gallbladder,  liver,  bone,  gill,  kidney,  and  ovary. 
Liver  neoplasia  was  not  a  significant  endpoint  in  any 
uninitiated  treatment  group. 

In  conclusion,  BDCM  tissue  levels  were  established  for 
short-term  exposed  fish  for  0,  1.5,  and  15  mg/L  BDCM.  A 
preliminary  PBPK  model  for  BDCM  in  the  medaka  was  developed. 
It  represents  the  first  PBPK  model  for  a  small  fish  that 
distinguishes  between  different  tissue  groups.  Data  to 
support  either  compartmental  or  PBPK  models  for  small  fish 
are  based  on  'lumped'  kinetic  data,  which  is  usually  the 
entire  fish.  The  experimental  design  of  this  project 
allowed  development  and  partial  validation  of  a  true  PBPK 
model  for  the  medaka.  This  modeling  effort  provided 
insights  into  important  species  extrapolation  issues  related 
to  exposure  and  dosimetry  of  BDCM. 

While  long-term  fish  exposure  to  two  concentration 
levels  of  BDCM  alone  did  not  cause  neoplasia,  significant 
findings  were  seen  at  the  cellular  level  in  the  gallbladder, 
liver,  bone,  gill,  kidney,  and  ovary  in  fish  dosed  with  15 
mg/L  BDCM.  Publication  of  this  work  in  the  open  literature 
is  pending. 
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BACKGROUND 

Previous  work  by  the  US  Army  Center  for  Environmental 
Health  Research  (USACEHR)  with  Japanese  medaka  ( Oryzias 
latipes)  exposed  to  0,  0.015,  0.15,  and  1.5  mg/L 
bromodichloromethane  (BDCM)  in  aquarium  water  for  nine 
months  showed  that  intrahepatic  concentrations  of  BDCM  were 
roughly  equivalent  to  the  concentration  of  BDCM  in  aquarium 
water  (Toussaint  et  al.,  2001c).  USACEHR  proposed  a  two- 
phase  study  to  determine  uptake  (Phase  I)  and  depuration 
(Phase  II)  kinetics  with  medaka  and  BDCM.  This  IAG  covered  • 
the  Phase  I  portion  of  the  project.  USACEHR  planned  to 
conduct  the  fish  exposures  and  perform  sample  analyses, 
while  an  extramural  collaborator  at  Wright  Patterson  Air 
Force  Base  built  the  physiologically  based  pharmacokinetic 
( PBPK)  model  for  medaka. 

PROJECT  DESCRIPTION 

PBPK  model  development:  Mature  Japanese  medaka  were 
statically  exposed  to  BDCM  <9  hours  at  three  BDCM 
concentrations.  The  proposed  number  of  fish  per  treatment 
level  was  20  males  and  20  females  (40  fish  per  treatment  x  4 
treatments  =  160  fish.  The  proposed  BDCM  treatment  levels 
were  0,  1.5,  15,  and  150  mg/L.  Blood  and  tissue  (gill, 
liver,  muscle,  spleen,  and  kidney)  samples  removed  from 
these  fish  were  to  be  chemically  analyzed  for  BDCM  and 
corrected  for  sample  volume.  All  raw  data  were  summarized 
in  spreadsheets  for  statistical  modeling. 

Six  month  bioassay:  Juvenile  Japanese  medaka  were 
exposed  to  two  concentrations  (1.5  and  15  mg/L)  of  BDCM  in  a 
continuous  flow-through  dosing  system.  Prior  to  BDCM 
dosing,  half  of  the  sample  groups  were  statically  exposed  to 
10  mg/L  diethylnitrosamine  (DEN)  for  48  hours.  This  level 
of  DEN  has  been  shown  to  alter  cells  (initiation)  but  does 
not  promote  neoplasia  (Hinton  et  al.,  1988).  After  28  days 
of  BDCM  dosing,  5  fish  were  removed  from  each  tank  for 
assessment  of  cell  proliferation.  The  remaining  25  fish  per 
replicate  continued  with  BDCM  dosing  throughout  the  six- 
month  test  duration.  At  test  termination,  the  fish  were 
grossly  observed  and  then  prepared  for  histopathology .  The 
pathology  results  were  analyzed  statistically  to  compare 
treated  fish  to  control  fish  at  the  tissue  level. 


INTRODUCTION 

BDCM  is  a  trihalome thane  that  is  frequently  found 
drinking  water  as  a  disinfection  by-product.  At  the  time 
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the  project  was  proposed,  BDCM  was  reported  as  a  mammalian 
liver  and  kidney  carcinogen  when  high  concentrations  (25  to 
50  mg/kg) were  administered  via  gavage  to  rodents (NTP,  1987). 
These  BDCM  concentrations  are  several  orders  of  magnitude 
higher  than  actual  drinking  water  levels  of  0.067  mg/L 
reported  by  Borum  et  al  (1998) .  A  need  existed  to  perform 
animal  studies  at  BDCM  concentrations  closer  to  drinking 
water  levels  (Boorman  et  al.,  1999). 

Alternative  animal  models  such  as  fish  are  a  natural 
choice  for  testing  aqueous  solutions.  Fish  have  been  shown 
to  be  sensitive  to  trace  levels  of  contaminants  in  aquatic 
media  (Gardner  et  al.,  1990;  Metcalfe  1989;  Twerdok  et  al . , 
1997;  Toussaint  et  al.,  2001b).  As  water-dwelling 
organisms,  fish  are  immersed  in  the  exposure  solution,  with 
dosing  occurring  through  dermal,  oral,  and,  most 
significantly,  through  respiratory  routes. 

Japanese  medaka  ( Oryzias  latipes)  are  hardy,  small  in 
size,  easy  to  culture,  and  have  a  relatively  short  time-to- 
tumor  response.  Medaka  are  reproductively  mature  at  three 
months  of  age,  and  an  exposure  duration  of  six  months 
encompasses  the  majority  of  the  active  growth  and 
development  phase  of  the  fish  life  cycle.  Sections  of  the 
entire  animal  will  fit  on  two  to  three  microscope  slides  so 
that  examination  of  each  tissue  is  possible  in  its 
anatomical  context.  The  low  rate  of  spontaneous  neoplasms 
in  medaka  (^0.3%)  aids  in  the  interpretation  of  bioassay 
results  (Hawkins  et  al.,  1995). 

Research  began  on  this  project  in  August  1999.  The 
animal  exposure  portion  of  the  experiment  at  USACEHR  ended 
in  May  2000,  while  statistical  modeling  was  completed 
December  2001.  Submission  of  this  work  for  publication  is 
planned  for  2002. 


MATERIALS  AND  METHODS 

Animal  Care.  Research  was  conducted  in  compliance  with 
the  Animal  Welfare  Act,  and  other  Federal  statutes  and 
regulations  relating  to  animals  and  experiments  involving 
animals  and  adheres  to  principles  stated  in  the  Guide  for 
the  Care  and  Use  of  Laboratory  Animals,  National  Academy 
Press,  Washington,  DC,  1996,  in  facilities  that  are  fully 
accredited  by  the  Association  for  the  Assessment  and 
Accreditation  of  Laboratory  Animal  Care,  International. 

Japanese  medaka  were  supplied  from  USACEHR  in-house 
cultures  and  were  reared  according  to  USACEHR  Standing 
Operating  Procedures,  as  described  previously  (Toussaint 
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2001a).  Fish  culture  conditions  were:  water  temperature 
25±1°C,  water  flow  of  100  mL/min  of  an  aerated,  hard, 
processed  groundwater,  an  18/6  hr  light/dark  cycle  with  full 
spectrum  fluorescent  lighting,  and  an  age-related  amount  of 
fish  flake  food  and  live  brine  shrimp  fed  daily. 

Test  Materials.  BDCM  (CAS  No.  75-27 -4)  of  99.8%  purity 
was  obtained  from  Aldrich  Chemical  Company,  Milwaukee,  WI. 
Under  test  conditions,  BDCM  did  not  break  down  to  form  other 
compounds,  as  demonstrated  by  chemical  analyses  of  aquaria 
water.  The  analysis  of  BDCM  for  aqueous  and  tissue  samples 
was  performed  using  a  flame  ionization  detector  interfaced 
to  a  model  6890  capillary  gas  chromatograph  (both  from 
Hewlett-Packard,  Wilmington,  DE)  equipped  with  a  model  7694 
automatic  static  headspace  sampler  and  fused  silica 
capillary  column  (30  m  x  0.25  mm  inner  diameter)  coated  with 
cross-linked  1%  methylsilicone  gum  phase,  film  thickness, 
0.33  urn  (Hewlett-Packard).  All  stocks  and  standard  solutions 
were  prepared  fresh  daily.  Sample  analysis  has  been 
previously  described  (Toussaint  et  al.,  2001a).  Triplicate 
5  mL  aliquots  of  the  40  mL  aqueous  samples  were  analyzed. 

The  detection  limit  for  this  gas  chromatography  method  was 
0.001  mg/L  and  the  average  percent  recovery  was  95%. 

DEN  (CAS  No.  55-18-5)  of  >99%  purity  was  obtained  from 
Sigma  Chemical  Company,  St.  Louis,  MO.  The  analysis  of  DEN 
was  performed  using  gas  chromatography  on  a  Hewlett  Packard 
(HP)  5880A  gas  chromatograph  with  a  flame-ionization 
detector  and  an  HP  19392A  integrator.  The  detection  limits 
were  <1  mg/L  and  the  percent  recovery  was  105%. 

BrdU  (CAS  No.  59-14-3)  of  >99%  purity  was  obtained  from 
Sigma  Chemical  Company,  St.  Louis,  MO.  The  analysis  of  BrdU 
was  performed  by  high  performance  chromatography  (HPLC)  on 
an  HP  1050  series  HPLC  equipped  with  variable  wavelength 
detector,  autosampler,  and  HP  3396A  integrator.  The 
detection  limit  was  1.77  mg/L  and  the  percent  recovery  was 
99.3%. 

Test  108-001:  Fish  Blood  and  Tissue  Analyses .  After  4 
to  9  hours  of  BDCM  exposure,  fish  were  euthanized  and 
tissues  harvested.  Fish  blood  was  collected,  liver,  gill, 
muscle,  kidney,  and  spleen  were  excised  from  each  animal, 
weighed,  sealed  in  10  mL  glass  headspace  vials,  and  then 
stored  at  0°C  until  analysis.  Due  to  their  small  size, 
spleens  and  kidneys  were  pooled  by  organ  type  within  each 
replicate.  Approximately  1500  tissue  samples  were 
chemically  analyzed  for  BDCM.  Thawed  samples  were 
individually  denatured  in  5  mL  of  2%  (w/v)  sodium  dodecyl 
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sulfate  at  60°C  for  2  hr.  Headspace  analysis  by  capillary 
gas  chromatography  immediately  followed  as  described  above. 
The  detection  limit  for  this  gas  chromatography  method  was 
0.001  mg/L  and  the  percent  recovery  was  95%. 


Test  108-001:  Test  Design.  For  Test  108-001,  the  age 
of  the  mature  fish  used  for  this  one-day  experiment  was  six 
months  of  age.  The  fish  were  not  fed  during  the  exposure. 
Prior  to  the  test,  the  population  of  fish  was  separated  into 
same  gender  groups.  Randomization  of  20  males  to  each  of  9 
test  vessels  and  of  20  females  to  each  of  another  9  test 
vessels  was  done  using  US ACE HR' s  computerized  "RANDOM" 
program.  A  summary  of  the  experimental  design  for  Test  108- 
001  is  shown  in  Table  I.  The  planned  duration  for  this 
experiment  was  4  hours . 


Table  I.  Test  108-001  Test  Design 


Group  ID 

Gender 

BDCM  Treatment  Level 

1 

male 

15  mg/L 

2 

male 

15  mg/L 

3 

male 

15  mg/L 

4 

female 

15  mg/L 

5 

female 

15  mg/L 

6 

female 

15  mg/L 

7 

male 

1.5  mg/L 

8 

male 

1.5  mg/L 

9 

male 

1.5  mg/L 

10 

female 

1.5  mg/L 

11 

female 

1 . 5  mg/L 

12 

female 

1.5  mg/L 

13 

male 

0  mg/L 

14 

male 

0  mg/L 

15 

male 

0  mg/L 

16 

female 

0  mg/L 

17 

female 

0  mg/L 

18 

female 

0  mg/L 

Test  108-001.  PBPK  Modeling  Methods.  Concentrations 
of  BDCM  were  modeled  in  various  tissues,  including  blood, 
liver,  muscle,  and  gill.  Remaining  tissue  was  lumped 
together  and  represented  as  'other'  compartment.  Metabolic 
parameters,  Vmax  and  Km,  for  metabolism  in  the  liver  were 
taken  from  the  rat  (Lilly  et  al.,  1998);  Vmax  was  scaled 
from  the  reference  value  in  the  rat  using  the  relationship 
of  reference  value  raised  to  the  0.75  power  and  multiplied 
by  body  weight  of  the  fish.  Cardiac  output  was  taken  from  a 
reference  value  in  a  one  kilogram  rainbow  trout  (Nichols  et 
al . ,  1996)  raised  to  the  0.75  power  and  multiplied  by  body 
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weight.  Blood  flow  through  tissues  was  determined  by 
multiplying  the  relative  volume  of  the  organ  by  the  body 
weight  scaled  cardiac  output;  thus,  blood  flow  through  the 
various  tissues  was  also  scaled.  Adjustments  to  the  model 
were  made  in  order  to  provide  better  fits  of  model-generated 
concentrations  of  BDCM  in  tissues  as  compared  to  assayed  or 
measured  BDCM  concentrations  in  tissues. 

Partition  coefficients  were  estimated  by  calculating 
the  concentration  ratio  of  tissue/blood  for  the  muscle, 
gill,  and  liver  in  the  low  dose  group  (1.5  mg/L  BDCM).  The 
water/blood  ratio  was  determined  using  measured  water 
concentrations.  The  'other'  compartment  represents  a 
lipophilic  compartment.  BDCM  has  modest  lipophilicity  in 
rat  fat  tissue  (Lilly  et  al.,  1998) .  In  the  medaka,  no  fat 
tissue  was  removed  for  analysis.  Our  modeling  efforts 
suggest  that  a  fat  depot  probably  exists  for  BDCM  in  the 
medaka.  We  were  unable  to  simulate  the  kinetic  behavior  of 
BDCM  in  the  medaka  fish  without  this  assumption.  Six  months 
of  effort  was  devoted  to  developing  experimental  methods  to 
measure  partition  coefficients  of  BDCM  in  small  pieces  of 
medaka  tissue  with  little  success. 

The  differential  equations  that  constituted  the 
structural  model  to  simulate  the  pharmacokinetics  of  BDCM 
are  shown  in  Figures  I  and  II.  A  full  list  of  all  the 
parameters  used  in  the  model  may  be  found  in  Table  II.  A 
good  fit  was  defined  as  model  predicted  BDCM  tissue 
concentrations  falling  within  one  standard  deviation  of  the 
measured  tissue  concentration.  The  model  simulation  needed 
to  be  robust  and  to  fit  all  tissue  concentration  profiles 
for  both  the  1.5  and  15  mg/L  BDCM  treatment  groups. 

Figure  I :  Differential  equations  used  to  model  the 
concentration  of  bromodichlorome thane  in  the  Japanese  medaka 
( Oryzias  latipes) 

dAG/dt  =  QC (CT  -  CA)  +  QP*GEE (CI*PGW  -  *CG/PGW) 

dAM/dt  =  QM (CA  -  AM/ (VM*PM) ) 

dAL/dt  =  QL (CA  -  (AL/ (VL*PL) )  -  RAMX) 

RAMX  =  dAMX  =  VMAX* (AL/ (VL*PL) ) / (KM  +  (AL/ (VL*PL) ) 
dAO/dt  =  QO (CA  -  (AO/VO*PO) ) 

CA  =  AG/VG*PB 

CV  =  (QM* (AM/VM*PM)  +  QL* (AL/VL*PL)  +  QO* (AO/VO*PO) ) /QC 
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Figure  II:  Structural  model  of  the  physiologically  based 
pharmacokinetic  model  for  bromodi chi orome thane  in  the 
Japanese  medaka  ( Oryzias  latipes) 


Model  terms  end  with  the  first  letter  of  the  tissue  with  which  they  are  associated  (G,  L,  M,  O,  A,  and  V  for 
gill,  liver,  muscle,  other,  arterial,  and  venous).  Model  terms  beginning  with  Q  are  flow  rates  (L/h)  such  that 
QC  is  cardiac  output,  QL  is  blood  flow  to  the  liver,  etc.  QP  is  rate  of  intake  into  the  gill.  Model  terms 
beginning  with  C  represent  concentrations  (mg/L);  arterial  concentrations  enter  tissues  and  leave  as  venous 
concentrations  (eg  CA  enters  liver  and  exits  as  venous  liver  concentration,  CVL).  Cl  and  CX  are  inspired 
and  expired  concentrations;  CT  is  a  transit  term  used  in  ACSL  code  and  is  equivalent  to  CV.  Partition 
coefficients  begin  with  P,  end  with  the  first  letter  of  the  tissue  with  which  they  are  associated. 


Test  108-002:  Test  Design.  For  Test  108-002,  the 
initial  age  of  the  medaka  was  16  days  +1  day,  and  the  fish 
were  fed  throughout  the  study.  Food  was  withheld  24  h  prior 
to  the  six-month  histopathology  sacrifice. 

Randomization  of  thirty  juvenile  fish  to  each  of  the  18 
test  vessels  was  done  using  USACEHR' s  computerized  "RANDOM" 
program.  A  summary  of  the  experimental  design  for  Test  108- 
002  is  shown  in  Table  III.  Under  the  continuous  flow¬ 
through  dosing  regimen,  sufficient  BDCM  remained  in  solution 
during  testing  to  maintain  dosage  levels.  Due  to  the  large 
number  of  test  aquaria  needed  for  the  test  design,  two 
solenoid-controlled  proportional  diluters  were  used  to 
deliver  test  material.  Stock  solutions  of  1785  and  1036 
mg/L  BDCM  were  made  daily  by  adding  neat  BDCM  to  processed 
well  water.  The  BDCM  stocks  were  stirred  for  24  +2  h  in 
sealed  glass  containers,  then  pumped  into  a  diluter  dosing 
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bottle.  Approximately  25  L  of  BDCM  stock  were  needed  for 
each  test  day. 

Table  II:  Model  parameters  and  description  for  the 
physiologically  based  pharmacokinetic-pharmacodynamic  model 
of  bromodichlorome thane  in  the  Japanese  medaka  (Oryzlas 
latipes) _ _ _ 


Scalable  rates 

QCC 

Cardiac  output 

2.07 

QPC 

Rate  of  uptake  at  gill 

7.06 

Fractional  blood  flows  (%CO) 

QLCb 

Fractional  blood  flow  to  liver 

0.03 

QMCb 

Fractional  blood  flow  to  muscle 

0.6 

QGC 

Fractional  blood  flow  to  gill 

i 

Relative  volumes  (%BW) 

VGC 

volume  of  gill 

0.1 

VLC  b 

volume  of  liver 

0.02 

VMC  b 

volume  of  muscle 

0.78 

Partition  Coefficients 

PGW 

Gill/water 

2.45 

PB 

Gill/blood 

0.46 

PL 

Liver/blood 

1.25 

PM 

Muscle  to  blood 

0.12 

POb 

“other’5  /  blood 

18 

Metabolism 

VMAXC  (mg/L/h) 

Scaleable  Vmax 

0.0128 

KM a  (mg/L) 

Michelis-Menten  constant  for 

metabolism 

0.5 

a  taken  (and  adjusted)  from  Lilly  et  al.,  1998 
b  Nichols  et  al.,  1996 


Six  test  aquaria  were  randomly  assigned  to  each  nominal 
concentration:  0,  1.5,  and  15  mg/L  BDCM.  The  test  began 
when  sixteen-day-old  fry  (+1  day)  were  randomized  to  each 
test  aquarium,  with  30  fry  per  five-gallon  aquarium.  The 
diluter  was  set  to  deliver  300  +15  mL  to  each  aquarium  every 
3  min  +15  seconds,  yielding  9-10  tank  volume  turnovers  per 
day. 

The  test  began  with  a  48-hour  exposure  to  0  or  10  mg/L 
DEN.  The  fish  were  rinsed  for  24  hours  in  processed  well 
water  before  beginning  the  six-month  BDCM  exposure.  After 
28  days  of  BDCM  exposure,  5  fish  per  aquarium  per  timepoint 
were  removed  for  hepatocellular  proliferation  assays.  The 
method  for  cell  proliferation  assessment  has  been  described 
previously  (Brennan  et  al.,  2001).  Briefly,  the  fish 
received  a  static  exposure  to  75  mg/L  5-bromo-2'- 
deoxyuridine  (BrdU)  for  72  hours,  followed  by  euthanasia, 
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fixation  in  acetone,  embedding,  sectioning,  staining,  and 
reading  of  the  slides.  The  remainder  of  the  test  fish 
received  six  months  of  BDCM  exposure,  followed  by  necropsy 
for  histopathological  evaluation. 


Table  III.  Test  108-002  Test  Design 


Group  ID 

BDCM  Treatment 

DEN  Treatment 

1 

0  mg/L 

0  mg/L 

2 

0  mg/L 

0  mg/L 

3 

0  mg/L 

10  mg/L 

4 

0  mg/L 

10  mg/L 

5 

1 . 5  mg/L 

0  mg/L 

6 

1.5  mg/L 

10  mg/L 

7 

15  mg/L 

0  mg/L 

8 

15  mg/L 

0  mg/L 

9 

15  mg/L 

10  mg/L 

10 

15  mg/L 

10  mg/L 

11 

0  mg/L 

0  mg/L 

12 

0  mg/L 

10  mg/L 

13 

1.5  mg/L 

0  mg/L 

14 

1.5  mg/L 

0  mg/L 

15 

1.5  mg/L 

10  mg/L 

16 

1.5  mg/L 

10  mg/L 

17 

15  mg/L 

0  mg/L 

18 

15  mg/L 

10  mg/L 

Test  108-002:  Chronic  Histopathology.  At  six  months 
of  BDCM  exposure,  fish  were  euthanized  with  an  overdose  of 
tricaine  methanesulfonate  (MS-222).  Fish  necropsy 
procedures,  fixation,  sectioning,  staining,  and 
histopathology  were  as  described  previously  (Toussaint  et 
al . ,  1999).  Briefly,  fish  were  opened  ventrally,  flushed 
with  Bouin' s  solution  and  sealed  in  tissue  cassettes.  After 
48  hr  immersion  in  Bouin' s,  the  tissue  was  rinsed  with  two 
successive  24  hr  treatments  of  70%  ethanol,  and  then  held  in 
10%  formalin.  Five  longitudinal  step  sections  were  cut  from 
paraffin  blocks,  stained  with  hematoxylin  and  eosin, 
followed  by  microscopic  evaluation  of  the  tissue.  The 
following  tissues  were  assessed  histologically:  bone 
(vertebra) ,  brain,  chromaffin  tissue,  corpuscle  of  Stannius, 
esophagus,  eye,  gallbladder,  gill,  heart,  hematopoietic 
tissue,  interrenal  tissue,  intestine,  kidney,  liver,  nares, 
ovary,  pancreas,  peripheral  nerve,  pineal  organ,  pituitary 
gland,  pseudobranch,  skeletal  muscle,  skin,  spinal  cord, 
spleen,  stato-acoustic  organ,  swim  bladder,  testis,  thymus, 
thyroid  tissue,  urinary  bladder,  and  gross  lesions. 

Test  108-002:  Statistical  Analyses.  Lengths,  weights, 
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and  histopathology  data  were  analyzed  statistically  (Snedcor 
and  Cochran,  1980;  Hsu  JC,  1996;  McCullagh  and  Nelder  1983) . 
Comparisons  of  length  and  weight  data  were  done  by  one-way 
analysis  of  variance  F  tests.  When  comparing  the  combined 
controls  versus  each  combined  dose-level  group  the  effect  of 
multiple  comparisons  was  accounted  for  by  Dunnet's  method. 
The  effect  of  multiple  comparisons  was  accounted  for  by 
Tukey' s  method.  For  each  replication  separately  and  for  the 
three  replications  combined,  standard  linear  regression 
methods  were  used.  For  histopathology  endpoints,  standard 
chi-square  tests  were  used  to  compare  the  occurrence  rates 
from  two  or  more  experimental  groups.  Fisher-exact  tests 
were  used  when  computationally  feasible.  The  effect  of 
multiple  comparisons  when  combining  two  or  more  experimental 
groups  was  accounted  for  by  the  Bonferroni  method.  For  each 
replication  separately  and  for  the  three  replications 
combined,  logistic  regression  methods  were  used  to  see  if 
the  occurrence  rate  of  each  endpoint  was  related  to  dose 
level . 

Archiving  Requirements.  Stained  fish  slides  were 
archived  and  will  be  held  for  a  minimum  of  five  years. 
Chemistry  samples  were  spent  during  analysis  and  were  not 
archived.  All  raw  data,  data  summaries,  and  test  reports 
will  be  maintained  by  USACEHR  for  a  minimum  of  10  years 
after  completion  of  the  test  report. 


RESULTS  AND  DISCUSSION 

Test  108-001:  Water  Quality.  Single  measurements  of 
temperature,  pH,  dissolved  oxygen,  and  conductivity  were 
taken  on  each  of  the  18  test  vessels  at  the  conclusion  of 
the  exposure.  Alkalinity  and  hardness  were  measured  in  one 
control  tank.  All  measurements  were  within  acceptable 
ranges.  Neither  pH  nor  dissolved  oxygen  levels  were  changed 
by  BDCM  concentration.  Results  were  25.1°C  for  temperature, 
7.3  for  pH,  7.5  mg/L  for  dissolved  oxygen,  645  umho/cm  for 
conductivity,  128  mg/L  for  alkalinity,  and  196  mg/L  for 
hardness . 

Test  108-001:  External  BDCM  Concentration.  For  Test 
108-001  (PBPK  test)  the  measured  BDCM  concentrations  were 
averaged  from  the  initial  and  final  readings  of  six 
replicate  containers  at  each  treatment  level  with  the 
following  mean  values  and  standard  deviations  obtained: 
below  detection  limits  of  0.001  mg/L  for  controls,  1.82  + 
0.2246  mg/L,  and  17.14  +  2.3552  mg/L. 

Test  108-001:  BDCM  in  Fish  Blood  and  Tissue.  At  the 
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time  of  fish  necropsy,  fish  were  removed  from  a  test 
container  and  then  held  in  clean  well  water  until 
euthanasia.  A  team  of  researchers  worked  on  each  container, 
with  an  average  working  time  of  approximately  70  minutes  per 
test  container.  Comparisons  of  measured  values  for 
individual  fish  tissues  within  each  group  of  the  beginning 
and  final  fish  necropsied  for  each  container  did  not  show  a 
consistent  depuration  trend,  and  it  is  assumed  that  this 
holding  period  had  minimal  impact  on  the  study. 

Additionally,  the  duration  of  the  exposure  had  a  wide 
variability  due  to  the  amount  of  time  required  for  necropsy. 
The  15  mg/L  fish  had  a  duration  range  of  6.5  to  8.5  hours, 
the  1.5  mg/L  fish  had  a  duration  range  of  5.0  to  8.5  hours, 
and  the  controls  had  a  duration  range  of  4.0  to  5.5  hours. 
Tissue  concentrations  for  the  1500  samples  are  shown  in 
Appendix  I,  with  raw  data  and  averages  by  group,  by  gender, 
and  by  treatment  level.  For  example,  female  medaka  weighing 
562  mg  +91  mg  (mean  and  standard  deviation)  were  found  to 
have  the  following  mean  BDCM  tissue  levels  at  15  mg/L  BDCM 
external  concentration:  8.77E+05  mg/L  in  blood,  2.15E+06 
mg/L  in  liver,  4.17E+05  mg/L  in  muscle,  1.70E+06  mg/L  in 
gill,  3.32E+06  mg/L  in  kidney,  and  4.45E+07  mg/L  in  spleen. 
Similar  order  of  magnitude  findings  occurred  in  males  at  the 
same  BDCM  treatment  level. 

Test  108-001:  PBPK  Model.  The  model  simulated 
concentrations  of  BDCM  in  medaka  fish  tissues  approached 
observed  concentrations  only  after  critical  adjustments  to 
the  PBPK  model  were  attempted  and  by  using  concentration 
ratios  from  high  dosed  male  medaka  as  partition 
coefficients.  The  model  code  is  given  in  Appendix  II,  while 
the  simulation  plots  are  shown  in  Appendix  III.  A  term  for 
gill  extraction  efficiency,  GEE,  was  added  to  the  PBPK  model 
in  the  gill  compartment  to  account  for  diffusional 
limitations  after  initial  model  runs  produced  highly 
exaggerated  tissue  concentrations.  The  value  for  GEE  was 
based  on  a  reported  extraction  efficiency  of  7%  for  butanol 
in  trout  gills  (McKim  et  al . ,  1995) .  This  was  one  of  the 
significant  findings  from  our  modeling  efforts.  Only  a  very 
small  fraction  of  the  'inhaled'  BDCM  dissolved  in  water  is 
actually  taken  up  into  systemic  circulation.  However,  it  is 
important  to  note  that  our  approach  to  describing  the 
gill/water  exchange  is  taken  from  approaches  used  for 
inhaled  gases  that  have  a  high  water  solubility  and  modest 
lipid  solubility  (Fisher  et  al.,  2000).  Other  fish  PBPK 
models  for  much  larger  fish  describe  the  gill/water  exchange 
by  more  rigorous  methods  (McKim  et  al.,  1995;  Erickson  and 
McKim,  1990) . 
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Initial  iterations  done  with  the  maximum  rate  of 
metabolism  set  at  the  value  reported  in  rats  (Lilly  et  al . , 
1998),  scaled  to  the  0.75  power  and  multiplied  by  body 
weight  of  the  fish  failed  to  simulate  the  observed  BDCM 
concentration  profiles  measured  in  the  medaka.  Scaling  the 
metabolic  rate  from  rats  to  the  medaka  fish  resulted  in 
extremely  rapid  clearance  of  DBCM  from  systemic  circulation. 
We  then  set  the  maximum  rate  of  metabolism,  Vmaxc,  to  zero 
and  increased  its  value,  in  an  iterative  fashion,  until 
suitable  fits  were  seen.  Adjustments  to  Vmaxc  were  carried 
out  merely  to  provide  better  fits  of  the  data.  To  determine 
the  metabolic  constants,  Vmax  and  Km,  for  oxidation  of  BDCM 
in  the  medaka  liver  additional  experimentation  are  required. 
Metabolism  of  BDCM  is  thought  to  be  predominantly  dependant 
on  CPY2E1  in  the  rat  (Allis  et  al . ,  2001).  However,  one 
report  suggests  that  this  isoform  is  not  expressed  in  the 
medaka  (Lipscomb  et  al.,  1997).  Other  isoforms  of  the 
cytochrome  p450  system  as  well  as  other  enzymes  may 
metabolize  BDCM  in  the  medaka. 

Hepatocellular  adenoma  incidence  rates  in  the  medaka 
fish  were  chosen  as  a  pharmacodynamic  endpoint  because  of 
reported  findings  of  BDCM  induced  liver  neoplasms  in  rats 
(NTP,  1987;  Dunnick  et  al.,  1987).  One  of  the  important 
'information  gaps'  in  the  use  of  medaka  as  an  alternative 
model  for  rodents  is  understanding  the  differences  in  routes 
of  exposure.  The  model  derived  dosimetry  parameter,  area- 
under-the-concentration-versus-time-curve  for  liver  (AUCL) , 
was  selected  for  regression  analysis  with  the  percent 
incidence  of  hepatocellular  adenoma  found  in  the  medaka 
study.  There  was  very  high  linear  correspondence  of  percent 
liver  incidence  with  AUCL  (r2=1.0).  Bioassay  studies  for 
the  medaka  were  carried  out  over  several  weeks,  while  the 
pharmacokinetic  study  was  a  4-hour  study.  The  PBPK  model 
was  developed  based  on  the  4-hour  study  and  used  to 
extrapolate  to  longer-term  exposures  that  were  used  in  the 
medaka  bioassay  study.  All  regression  parameters  relating 
percent  incidence  to  the  exposure  metric  of  AUCL  are 
presented  in  Table  IV.  More  pharmacokinetic  data  is  needed 
to  describe  the  pharmacodynamic  responses  with  a  Hill 
equation  or  Weibull  cumulative  density  function,  which  would 
provide  more  a  more  accurate  description  of  the  dose-toxic- 
effect  relationship  of  BDCM  (Gibaldi  and  Perrier,  1982; 
Newman,  1995) . 
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Table  IV:  Pharmacodynamic  regression  parameters  -  percent 


effect  as  a  function  of  AUCL 


Exposure  metric 

slope 

intercept 

AUCL 

1  x  10-6 

1.48 

1.00 

Weekly  AUCL 

4  x  10-5 

1.48 

1.00 

Daily  AUCL 

2  x  10-4 

1.48 

1.00 

With  further  research,  improvements  in  the  PBPK-PD 
model  can  be  accomplished  in  several  areas.  With  more 
pharmacokinetic  studies,  a  richer  time  course  data  set  could 
be  used  to  obtain  optimized  model  parameters  using  Advanced 
Continuous  Simulation  Language  (ACSL)  optimization  routine. 

The  physiological  parameters  we  used  for  the  medaka  were 
scaled  down  from  large  fish.  Further  research  on  the 
physiology  (blood  flows  and  breathing  rates)  of  the  medaka 
would  improve  the  PBPK  model  robustness.  Separate  metabolic 
studies  are  required  to  quantify  the  metabolic  capability  of 
the  medaka.  Finally,  research  to  quantify  the  diffusional 
limitations  associated  with  gill  mucous  and  gill  diffusion 
would  allow  for  more  accurate  estimates  of  chemical  intake. 

Test  108-002:  Test  Water  Quality.  Water  quality 
parameters  were  measured  in  each  test  aquarium  on  the 
following  frequency:  daily  (temperature) ,  weekly  (dissolved 
oxygen,  pH,  conductivity,  alkalinity,  and  hardness) ,  and 
monthly  (ammonia) .  All  parameters  were  found  to  be  within 
acceptable  ranges.  Summarized  mean  values  with  standard 
deviations  are  24.8  +  0.4°C  for  temperature,  7.4  +  3.2  mg/L 
for  dissolved  oxygen,  666  +  108  umho/cm  for  conductivity, 
0.001  +  0  mg/L  for  un-ionized  ammonia.  The  pH  range  was 
from  7.2  to  8.0  standard  units. 

Test  108-002:  Cell  Proliferation.  Summarized  raw  data 
from  both  count  labeling  indices  (CLI)  and  area  labeling 
indices  (ALI)  are  shown  in  Table  V.  Cell  proliferation  was 
not  a  significant  finding  when  BDCM-treated  fish  were 
compared  to  control  fish.  For  both  ALI  and  CLI,  there 
appeared  to  be  a  linear  effect  of  BDCM  and  the  slope 
appeared  to  be  the  same  for  both  DEN=0  mg/L  and  DEN=10  mg/L 
and  there  appeared  to  be  no  effect  of  DEN  adjusted  for  BDCM. 
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Table  V.  Area  and  Count  Labeling  Indexes  Statistics  by  Tank 


Tank 

DEN 

BDCM 

N 

area  labeling 
mean  std 

count 

mean 

labeling 

std 

1 

0 

0 

5 

0.0022 

dev 

0.0032 

0.0772 

dev 

0.1094 

2 

0 

0 

5 

0.0032 

0.0034 

0.0856 

0.0867 

3 

10 

0 

5 

0.0266 

0.0247 

0.5006 

0.4436 

4 

10 

0 

5 

0.0270 

0.0315 

0.5382 

0.6415 

5 

0 

1.5 

5 

0.0150 

0.0110 

0.2240 

0.1367 

6 

10 

1.5 

5 

0.0098 

0.0086 

0.2268 

0.1844 

7 

0 

15 

5 

0.0244 

0.0388 

0.4634 

0.7192 

8 

0 

15 

5 

0.0788 

0.1183 

1.2356 

1.3473 

9 

10 

15 

5 

0.1502 

0.2193 

2.0978 

2.6063 

10 

10 

15 

5 

0.0392 

0.0528 

0.3766 

0.5084 

11 

0 

0 

5 

0.0112 

0.0116 

0.2490 

0.1888 

12 

10 

0 

5 

0.0048 

0.0054 

0.1134 

0.1101 

13 

0 

1.5 

5 

0.0134 

0.0128 

0.2066 

0.1655 

14 

0 

1.5 

5 

0.0254 

0.0241 

0.5386 

0.4974 

15 

10 

1.5 

5 

0.0204 

0.0206 

0.2974 

0.2671 

16 

10 

1.5 

5 

0.0224 

0.0257 

0.4412 

0.5684 

17 

0 

15 

5 

0.0226 

0.0318 

0.5542 

0.6573 

18 

10 

15 

5 

0.0268 

0.0220 

0.7430 

0.6261 

Test  108-002:  BDCM  Chemical  Analyses .  Over  the  course 
of  six  months,  weekly  chemical  analyses  were  performed  on 
each  aquarium.  Mean  measured  concentrations  (replicates 
combined)  were  obtained  from  averaging  the  27  weekly  BDCM 
measurements  of  all  six  test  aquaria  per  treatment  level  to 
obtain  the  following  mean  values  and  standard  deviations: 
0.008  +  0.017  mg/L  (controls),  1.468  +  0.379  mg/L,  and 
14.885  +  3.178  mg/L. 

Test  108-002:  Fish  Growth  and  Survival.  A  summary  of 
the  fish  growth  measurements  is  shown  in  Table  VI. 

Mortality  through  the  six-month  exposure  period  was 
acceptable  (<4%)  in  all  control  and  treated  groups.  A  one¬ 
way  analysis  of  variance  was  used  to  see  if  each  of  the  3 
groups  of  fish  at  each  treatment  level  in  each  group  had  the 
same  mean  length  and  weight  values.  In  general  this 
appeared  to  be  true  except  for  the  mean  length  among  the 
DEN=0/BDCM=0  and  DEN=0/BDCM=1 . 5  groups.  For  weight,  within 
each  treatment  group,  the  3  groups  of  fish  were  not 
significantly  different  from  each  other;  for  length, 
however,  equality  of  the  3  tanks  (25.2,  24.0  and  23.6)  of 
the  control  group  (DEN=0  BDCM=0)  and  the  3  tanks  (26.2,  24.8 
and  24.6)  of  the  treatment  group  DEN=0  and  BDCM=1 . 5  were 
significantly  different  from  each  other.  Because  the 
differences  in  length  for  the  3  replicate  control  groups  and 
for  the  3  replicates  of  the  treatment  group  DEN=0/BDCM=1 . 5 
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were  not  substantial,  the  3  replicate  groups  of  fish  were 
combined.  For  length,  there  appeared  to  be  a  BDCM  effect  for 
both  DEN  groups,  however  only  the  DEN=10  group  showed  a 
clear  linear  trend  (length  decreases  as  BDCM  increases) ; 
there  was  a  consistent  DEN  effect  (DEN=10  is  longer)  but  the 
statistical  significance  is  strong  for  the  BDCM=0  group, 
moderate  for  the  BDCM=1 . 5  group,  and  nil  for  the  BDCM=15 
group.  For  weight,  there  appeared  to  be  a  linear  BDCM 
effect  only  in  the  DEN=10  group  (weight  decreased  while  BDCM 
increased) ;  there  appeared  to  be  a  consistent  DEN  effect 
(DEN=10  is  heavier  than  DEN=0) ,  however,  the  statistical 
significance  of  this  effect  was  strong  in  the  BDCM=0  group 
while  nil  in  the  BDCM=1 . 5  and  BDCM=15  groups. 


Table  VI.  Test  108-002  Fish  Growth  Measurements 


BDCM 

DEN 

Mean  Length 

Mean  Weight 

Treatment 

Treatment 

(mm) 

(mg) 

Controls 

0  mg/L 

24.3 

280 

Controls 

10  mg/L 

25.5 

303 

1.468  mg/L 

0  mg/L 

25.2 

289 

1.468  mg/L 

10  mg/L 

25.6 

292 

14.885  mg/L 

0  mg/L 

23.3 

273 

14.885  mg/L 

10  mg/L 

23.7 

279 

Test  108-002:  Fish  Histopathology.  At  six  months  of 
BDCM  exposure,  ANOVA  indicated  that  the  significant  findings 
(p  =  0.05)  were  present  in  multiple  organ  systems,  as  shown 
in  Appendix  IV.  DEN  initiation  increased  findings  of 
hepatocellular  adenoma  were  only  observed  in  15  mg/L  males. 
Hepatocellular  carcinoma  was  not  significant  for  either  DEN 
treatment  or  BDCM  treatment,  regardless  of  DEN  treatment. 

The  finding  of  liver  basophilic  foci  was  significant  for  15 
mg/L  BDCM-exposed  fish  for  both  males  and  females.  Reports 
in  the  literature  indicate  that  basophilic  foci  in  medaka 
are  often  likely  to  develop  into  hepatocellular  adenomas 
(Bunton  1996;  Hinton  et  al.,  1988).  For  the  present  study, 
this  was  confirmed  with  15  mg/L  males  but  not  for  females. 
Liver  neoplastic  findings  were  not  significant  in  either 
BDCM  treatment  level  in  uninitiated  fish.  These  results 
suggest  that  BDCM  does  not  act  as  a  cancer  initiator  at 
either  treatment  level. 

Findings  of  bile  duct  concretions,  dilatation,  and 
hyperplasia  observed  in  a  previous  USACEHR  medaka  exposure 
to  BDCM  at  1.4  mg/L  (Toussaint  et  al . ,  2001c)  were  confirmed 
in  the  current  study,  while  gallbladder  concretions  were  not 
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significant.  Additionally,  gill  lamellae  fusion  was 
significant  for  all  1.4  mg/L  BDCM-exposed  medaka  except 
males  which  had  not  been  initiated  with  DEN.  Males  without 
DEN  initiation  showed  significant  levels  of  kidney  tubular 
casts . 

When  BDCM  treatment  was  increased  to  15  mg/L,  more 
statistically  significant  histology  endpoints  were  seen. 
Gallbladder  concretions  were  significant  for  both  males  and 
females,  while  gallbladder  epithelium  hyperplasia  was 
significant  for  uninitiated  females  and  for  DEN-initiated 
males.  Within  the  liver,  bile  duct  concretions  were 
significant  for  both  males  and  females  regardless  of  DEN 
treatment,  while  bile  duct  dilatation  and  bile  duct 
hyperplasia  were  significant  for  most  test  groups.  Bone 
malformations  were  common  in  15  mg/L  BDCM-treated  fish,  with 
the  jaw  and  vertebrae  significantly  affected.  Findings  in 
the  gill  included  fusion  of  gill  lamellae,  hyperplasia  of 
gill  epithelium,  malformed  arches,  and  malformed  filaments. 
Significant  levels  of  thyroid  follicles  were  only  seen  in 
DEN-initiated  males.  Within  the  kidney,  significant 
findings  of  tubular  casts,  tubular  degeneration,  and 
mineralization  occurred  independent  of  DEN  treatment,  but 
varied  by  gender.  Lastly,  ovarian  findings  of  early 
vitellogenic  stage  and  undeveloped  ovaries  were  significant 
for  females  regardless  of  DEN  treatment. 

Concretion  occurrence  in  gallbladder  and  bile  ducts  has 
previously  been  identified  in  BDCM-treated  medaka  (Toussaint 
et  al.,  2001c).  Gallstones  of  unknown  significance  or 
etiology  have  also  been  identified  in  tilapia  ( Oreochromic 
mossambicus)  (Ferguson,  1989) .  The  occurrence  of  biliary 
concretions  has  been  noted  in  monkeys,  cattle,  and  pigs, 
with  an  origin  of  a  crystallization  nidus  (Jones  and  Hunt, 
1983) .  It  was  suggested  that  a  nidus  might  result  from  an 
infectious  agent  such  as  Pseudomonas  aeruginosa ,  which  has 
been  identified  as  the  causative  agent  of  concretions  in 
sheep,  or  from  particles  of  bile  salts,  calcium  carbonate, 
cholesterol,  or  precipitated  bilirubin.  The  significance  of 
these  endpoints  and  other  chronic  pathology  findings  in 
BDCM-treated  medaka  is  unknown. 

During  the  last  six  weeks  of  the  chronic  BDCM  exposure, 
fish  behavior  in  the  15  mg/L  tanks  became  noticeably 
different  from  normal  behavior.  Observations  of  behavior  by 
fish  aquaria  were  noted  and  recorded  in  approximately  one- 
week  intervals.  The  control  fish  showed  normal  swimming  and 
feeding  behavior.  They  were  distributed  evenly  throughout 
the  water  column.  The  1.5  mg/L  BDCM  fish  showed  normal 
feeding  behavior,  although  they  seemed  to  have  a  heightened 
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response  to  stimuli.  These  fish  fed  throughout  the  water 
column.  At  15  mg/L  BDCM,  abnormal  fish  behaviors  were 
observed.  All  fish  had  rapid  respiration  and  were 
hyperexcitable .  When  movement  was  observed,  it  was 
exaggerated  and  without  regard  to  fish  orientation,  i.e., 
the  fish  would  veer  off  on  a  90°  angle  towards  the  top  or 
bottom  of  the  aquaria.  These  fish  were  not  evenly 
distributed  throughout  the  water  column  and  were  seen 
hovering  at  the  top  or  resting  on  the  bottom.  Feeding  only 
occurred  at  the  water  surface.  Severe  abdominal  edema  was 
seen  in  two  fish.  Spinal  abnormalities  such  as  scoliosis  and 
kyphosis  were  observed,  along  with  jaw  abnormalities.  DEN 
treated  fish  exhibited  the  same  behaviors  as  untreated  fish 
at  each  BDCM  treatment  level.  Since  these  observations  were 
not  blinded,  the  data  were  not  statistically  analyzed; 
however  they  confirm  that  the  15  mg/L  fish  were  most 
effected. 


CONCLUSIONS 

BDCM  tissue  levels  were  established  for  short-term 
exposed  fish  for  0,  1.5  and  15  mg/L.  A  preliminary  PBPK 
model  in  medaka  was  developed.  This  PBPK  modeling  effort 
represents  the  first  PBPK  model  for  a  small  fish  that 
distinguished  between  different  tissue  groups  and  provides 
insights  into  species  extrapolation  issues  related  to 
exposure  and  dosimetry  of  BDCM.  Chronic  exposure  to  1.5  or 
15  mg/L  BDCM  caused  significant  changes  in  fish  tissue 
pathology  after  six  months,  with  gallbladder,  liver,  bile 
ducts,  kidney,  ovary,  gill,  and  bone  having  pathological 
effects.  DEN  initiation  increased  the  occurrence  of 
hepatocellular  adenomas  in  males  and  basophilic  foci  in  both 
genders  at  the  15  mg/L  BDCM  treatment  level. 
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Appendix  I.  BDCM  Tissue  Levels  in  Medaka  (shaded  areas  were  below  detection  limits  and  changed  to  the  0.001  detection  limit  for  calculation  purposes) 
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Appendix !.  BDCM  Tissue  Levels  in  Medaka  (shaded  areas  were  below  detection  limits  and  changed  to  the  0.001  detection  limit  for  calculation  purposes) 
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Appendix  I.  BDCM  Tissue  Levels  in  Medaka  (shaded  areas  were  below  detection  limits  and  changed  to  the  0.001  detection  limit  for  calculation  purposes) 
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Appendix  I.  BDCM  Tissue  Levels  in  Medaka  (shaded  areas  were  below  detection  limits  and  changed  to  the  0.001  detection  limit  for  calculation  purposes) 
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Appendix  I.  BDCM  Tissue  Levels  in  Medaka  (shaded  areas  were  below  detection  limits  and  changed  to  the  0.001  detection  limit  for  calculation  purposes) 
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Appendix  I.  BDCM  Tissue  Levels  in  Medaka  (shaded  areas  were  below  detection  limits  and  changed  to  the  0.001  detection  limit  for  calculation  purposes) 
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Appendix  II:  Model  code 


Program:  medpdl  in  Medaka 
!  Jason  L  Boyd  and  J  W  Fisher,  PhD 
!  13  Nov  2001 

(metabolism  -  MM  in  liver,  based  on  rat  values  (Nichols  etal,  1998) 

!(Km(‘scaling’  Km  effectively  decrease  affinity),  scaled  Vmax) 

(partitions  -  using  "concentration  ratios"  from  low  dose  medaka  males 

(blood  flow  -  gill  compartment  sans  segmental  (parallel)  gill  blood  flow,  only  serial  blood  flow 

!gill  volume  QGC=1.0 

(subtracting  fractional  gill  volume  (0. 1)  from  !the  muscle  tissue  compartment  which  Nichols  etal  1993  call 
!  Viscera' 

!GEE  is  set  to  0.07 
!QPC  7.2,  Nichols  etal,  1996 

!  [Nichols  etal,  1996  compares  channel  cats  &  rainbow  trout  PBTK] 

!  -not  presented  here,  may  provide  better  description  of  flux  at  gill 
!~see  McKim  etal  1985 

j - initial  section  for  constants - 

INITIAL 

CONSTANT  QPC=7.06  Irate  of  uptake  at  the  gill  (L/(kgh))(Nichols  etall993,  Qw) 

CONSTANT  QCC=2.07  Icardiac  output  (L/(kgh))(Nichols  etal  1991, 1993) 

CONSTANT  QL00.03  Ifraction  blood  flow  to  liver  (Nichols  etal  1991,1993) 

CONSTANT  QMC=0.6  Ifraction  blood  flow  to  muscle  (Nichols  etal  1991, 1993) 

CONSTANT  QGC=1.0  Ifraction  blood  flow  to  gills 

CONSTANT  PGW=24  (partition  gill/water 

CONSTANT  PB=4.62  (partition  gill/blood 

CONSTANT  PL=1.3  I  partition  liver/blood 

CONSTANT  PM=0.12  Ipartition  muscle/blood 

CONSTANT  PO=  1 8  Ipartition  other/blood  (from  rat  fat) 

CONSTANT  BW=0.0005  Ibw  of  individual  fish  (kg) 

CONSTANT  VLC=0.015  (fractional  volume  of  liver  tissue  (Nichols  etal  1991/3) 

CONSTANT  VMC=0.781  (fractional  volume  of  muscle  tissue  (Nichols  etal  1991/3) 

CONSTANT  V000.094  (fractional  volume  of  'other*  tissue  (subtraction) 

CONSTANT  VGOO.  10  (fractional  volume  of  gill  tissue  (estimation) 

CONSTANT  VMAX012.8  (max  rate  of  hepatic  metabolism  (mg/hkg) 

CONSTANT  KMC=0. 5  (Michelis-Menten  rate  constant  (mg/L) 

CONSTANT  KFC=0.000000000001  !  1st  order  rate  constant  (h-1) 

CONSTANT  CONC=1.5  (cone  in  water  (mg/L) 

CONSTANT  GEE=0.07  (gill  extraction  efficiency  water  to  gill 

'Timing  commands’ 

CONSTANT  CINT=0.05 
CONSTANT  TSTOP=8.5 

CONSTANT  TCHNG=4  (time  (h)  of  exposure 

'SCALED  PARAMETERS' 

QC=QCC*BW**0.75  !(L/h)  bwscaled  cardiac  output 

QP=QPC*BW**0.75  !  (L/h)  bwscaled  rate  of  uptake  at  gill 

QL=QLC*QC 

QM=QMC*QC 

QG=QGC*QC 

QO=QC-QL-QM  !QO  the  rate  of  bl  flow  left  over  to  go  to  ’other'  tissues 
VL=VLC*BW 
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VM=VMC*BW 

VO=VOC*BW 

VG=VGC*BW 

KF=KFC*BW**0.75  INichols  etall993  assumed  no  biotransf 
VMAX=VMAXC*BW**0.75  INichols  etal,  1998 
KM=KMC  INichols  etal,  1998 

END  I 'END  OF  INITIAL' 

! - -DYNAMIC  SECTION- - - - 


DYNAMIC 

ALGORITHM  IALG=2  12  IS  GEAR'S  STIFF 

! - -DERIVATIVE  SECTION- - 

DERIVATIVE 

I - -'INHALATION'  SECTION- - 


procedural 

pflag  =  pulse(0.,20.,tchng)  llnhalation  switch  for  on  or  off 

Cl  =  CONC*pflag  llnhaled  cone 

end 

j - model  code - 

'AG  =  amount  in  gill  (mg)' 

RAG  =  QC*(CT  -  CA)  +  QP*(CI*PGW*GEE  -  CX) 

IRAG  =  QC*(CT  -  CA)  +  QP*((CI  -  CX)/((RM/PGW)  -  RAQ) 

I  SEE  YALKOWSKI  &  MOROZOWICH  1980  AND 
IMCKIM  ETAL  1985 
AG  =  INTEG(RAG,0.) 

CA  =  AG/(VG*PB) 

CG  =  AG/VG 
AUCG=  INTEG(CG,0.) 

AUCA=  INTEG(CA,0.) 

RAJ  =  QP*CI*PGW  I  this  is  rate  of  amount  input  into  gill  from  water 

AI  =  INTEG(RAI,0.) 

CX  =  GEE+CG/PGW 
RAX  =  QP*CX 
AX  =  INTEG(RAX,0.) 

DOSEINH  =  AI-AX  I  a  measure  of  extraction  efficiency 


'AM  =  amount  in  muscle  (mg)' 
RAM  =  QM*(CA-CVM) 

AM  =  INTEG(RAM,0.) 

CVM  =  AM/(VM*PM) 

CM  =  AM/VM 
AUCM=  INTEG(CM,0.) 


'AL  =  amount  in  liver  (mg)' 
RAL  =  QL*(CA-CVL)-RAMX 
AL1  =  INTEG(RAL,0.) 

AL  =  MAX(AL1,0.) 

CVL  =  AL/(VL*PL) 

CL  =  AL/VL 
AUCL  =  INTEG(CL,0.) 
AUCLW  =  AUCL/168 
AUCLM  =  AUCL/720 
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!the  AL  and  AL1  statement  is  an  attempt  to  fix  the  concentration  from  going  below  zero 

'AMX  =  amount  metabolized  (mg)' 

RAMX  =  VMAX*CL/(KM+CL) 

AMX  =  INTEG(RAMX,0.) 

AMXW  =  AMX/168 
AMXM  =  AMX/720 
CMX  =  AMX/VL 
AUCX=  INTEG(CMX,0.) 

AUCXW  =  AUCX/168 
AUCLM  =  AUCX/720 

'AO  =  amount  in  other  tissue  compartment  (mg)' 

RAO  =  QO*(CA-CVO) 

AO  =  INTEG(RAO,0.) 

CVO  =  AO/(VO*PO) 

CO  =  AO/VO 
AUCO  =  INTEG(CO,0.) 

'CV  =  mixed  venous  blood  concentration  (mg/L)' 

CV  =  (QM*CVM  +  QL*CVL  +  QO*CVO)/QC 
CVAUC=  INTEG(CV,0.) 

'CT  =  concentration  (mg/L)  while  in  transit  arterially  from  heart  to  gill' 

CT  =  CV 

'TMASS  =  mass  balance  (mg)' 

TMASS  =  AL+AM+AO+AX+AMX 

'DOSEX  =  net  amount  absorbed  (mg)' 

DOSEX  =  AI-AX-AMX  !I  added  AMX  to  this  because  AMX  would  be  leaving  also 

TERMT(T.  GE.TSTOP) 


END  !  'end  of  derivative' 

END  ! 'end  of  dynamic' 

END  !  'end  of  program' 


!  File:  med.cmd 

!  medaka  pbpk  cmd 

!  Jason  L.  Boyd  and  Jeff  Fisher,  PhD 

!  procedures:  3  sets  of  different  partitions,  automatically  run  and  plot 
procedure  parthm 

s  pgw=2.45,  pb=2.31,  pl=3.96,  pm=0.42 
s  title(l)='BDCM  in  CV,  CL,  high  cone  ratios' 
end 

Ido  not  put  a  semicolon  in  title,  semicolon  truncates  title 
procedure  partlm 

s  title(l)='BDCM  in  CV,  CL,  low  cone  ratios' 
s  pgw=2.45,  pb=0.46,  pl=1.25,  pm=0.12 
end 
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Ifirst  line  mean,  2nd  line  -sd,  3rd  line  +sd 
DATA  mhi  & 

(t,  cv,  cl,  cm,  eg) 

4.01,  15.9,  63.0, 6.66,  36.7 
4.01,  5.8,  0.0001,  2.43,  0.0001 
4.01,  26.0, 169,  10.9, 74.9 
end 

Ifirst  line  mean,  2nd  line  -sd,  3rd  line  +sd 
DATA  mlo  & 

(t,  cv,  cl,  cm,  eg) 

4.01, 7.97,  10.0, 0.96,  3.68 
4.01,  3.32,  3.79, 0.348,  1.8 
4.01,  12.6,  16.2, 1.57,  5.56 
end 

Ifollowing  data  excluded  in  preference  to  total  averaged  lump  data 
Idata  for  males  high  dose 
IDATA  mhigh  & 

!(t,  cv,  cl,  cm,  eg) 

14.017, 0.67,  2.54,  0.53, 2.7 
14.25, 0.60, 1.82,  0.28, 1.42 
15.0,  0.535, 1.46, 0.19,  0.99 
lend 
! 

Idata  for  males  low  dose 
IDATA  mlow  & 

!(t,  cv,  cl,  cm,  eg) 

14.017, 0.31,  0.22,  0.06, 0.18 
14.25, 0.29,  0.23,  0.06, 0.18 
15.0,0.29,0.23,0.04,0.14 
lend 
! 

Idata  for  females  high  dose 
IDATA  fhigh& 

!(t,  cv,  cl,  cm,  eg) 

14.017,1.07,3.67,  0.92,4.99 
14.25, 0.77,  2.80,  0.49, 1.87 
15.0, 0.832,  1.13, 0.30,  0.92 
lend 
! 

Idata  for  females  low  dose 
IDATA  flow  & 

!(t,  cv,  cl,  cm,  eg) 

14.017,0.43,0.69,0.14,0.41 
14.25, 0.42,  0.44,  0.14,  0.44 
15.0,  0.41, 0.42, 0.12,  0.35 
lend 

(previous  'data'  block  excluded  in  favor  of  analyzing  data  as  end  of  exposure,  lumped  data 

s  weditg=.false.  I  suppress  output  from  schedule 
s  hvdpm=  .false.  Isuppress  high  volume  display 

prepare  t,cv, cl, cm, eg, co,cvauc,aucl,dosex 

procedure  pith 
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start 

plot  /data=mhi  cv  /lo=-100  /hi=200  /tag='(mg/L)'  /char=2,  cl  /char=3,  cm  /chai=4,  eg  /char=5  /same 
end 

procedure  pltl 
start 

plot  /data=mlo  cv  /lo=-10  /hi=50  /tag='(mg/L)'  /char=2,  cl  /char=3,  cm  /char=4,  eg  /char=5  /same 
end 

procedure  hi  run 

s  conc=15,  title(37)='conc=15.0  mg/L' 
parthm 

s  vmaxc=0,  title(52)='vmaxc=0.0000' 
pith 

s  vmaxc=0.128,  title(52)='vmaxc=0.0128' 
pith 

s  vmaxc=0.128,  title(52)='vmaxc=0. 1280' 
pith 

s  vmaxc=1.28,  title(52)='vmaxc=1.2800' 

pith 

partlm 

s  vmaxc=0,  title(5  2)=' vmaxc=0 .0000' 
pith 

s  vmaxc=0.128,  title(52)='vmaxc=0.0128' 
pith 

s  vmaxc=0. 128,  title(52)='vmaxc=0. 1280' 
pith 

s  vmaxc=1.28,  title(52)=Vmaxc=1.2800' 

pith 

end 

procedure  lorun 

s  conc=1.5,  title(37)='conc=1.50  mg/L1 
parthm 

s  vmaxc=0,  title(52)='vmaxc=0.000' 

pltl 

s  vmaxc=0. 128,  title(52)='vmaxc=0. 128' 
pltl 

s  vmaxc=1.28,  title(52)='vmaxc=1.280' 
pltl 

s  vmaxc=12.8,  title(52)=' vmaxc=12.80' 

pltl 

partlm 

s  vmaxc=0,  title(52)='vmaxc=0.000' 
pltl 

s  vmaxc=0.128,  title(52)=‘vmaxc=0.128' 
pltl 

s  vmaxc=1.28,  title(52)='vmaxc=1.280' 
pltl 

s  vmaxc=12.8,  title(52)=Vmaxc=12.80' 

pltl 

end 


procedure  pol8 
s  po=18,  title(65)='po=18' 
hirun 
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lorun 

end 

procedure  po25 

s  po=25,  title(65)- po=25' 

hirun 

lorun 

end 

procedure  po28 

s  po=28,  title(65)- po=28f 

hirun 

lorun 

end 

procedure  po30 

s  po=30,  title(65)- po=30' 

hirun 

lorun 

end 

procedure  po23 

s  po=23,  title(65)-po=23' 

hirun 

lorun 

end 

procedure  po21 

s  po=21,  title(65)- po=21’ 

hirun 

lorun 

end 

procedure  run2 

s  conc=15,  title(37)-  conc=15.0  mg/L' 
s  vmaxc=0,  title(52)=,vmaxc=0.0000' 
pith 

s  conc=1.5,  title(37)=,conc,1.50  mg/L' 
s  vmaxc=0,  title(52)=,vInaxc=0.0000, 
pltl 

s  conc=15,  title(37)- conc=15.0  mg/L' 
s  vmaxc=0.0128,  title(52)=Vmaxc=0.0128, 
pith 

s  conc=1.5,  title(37)='conc=L50  mg/L‘ 
s  vmaxc=0.0128,  title(52)=,vInaxc=0.0128, 
pltl 

s  conc=15,  title(37)='conc=15.0  mg/L’ 
s  vmaxc=0.128,  title(52)=,vmaxc=0.1280' 
pith 

s  conc=1.5,  title(37)=,conc'1.50  mg/L’ 
s  vmaxc=0. 128,  title (5 2)= Vmaxc=0 . 1280' 
pltl 

s  conc=15,  title(37)=,conc=15.0  mg/L' 
s  vmaxc=1.28,  title(52)- vmaxc=L2800, 
pith 

s  conc=1.5,  title(37)='conc'L50  mg/L' 
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s  vmaxc=1.28,  title(52)=,vmaxc=1.2800' 
pltl 

s  conc=15,  title(37)='conc=15.0  mg/L' 
s  vmaxc=12.8,  title(52)=* vmaxc=12.800' 
pith 

s  conc=1.5,  title(37)='conc'1.50  mg/L' 
s  vmaxc=12.8,  title(52)=,vmaxc=12.800' 
pltl 
end 

procedure  bycoef 

parthm 

run2 

partlm 

run2 

end 

procedure  bypo 
s  po=18,  title(65)='po=18' 
bycoef. 

s  po=21,  title(65)='po=21' 
bycoef 

s  po=25,  title(65)='po=25' 
bycoef 

s  po=28,  title(65)='po=28' 

bycoef 

end 

procedure  doit 

pol8 

po30 

end 
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I:  Simulation  at  low  concentration  using  high  dosed  male  concentration  ratios  as  partitions  (same  as  plot  I,  except 
includes  the  “other”  compartment) 


Appendix  III.  PBPK  Plot  Simulations 


V:  High  dosed  male  medaka  using  high  dose  male  concentration  ratios  as  partition  coefficients,  and  Vmaxc 
0.0128.  (same  as  Plot  III  but  without  the  “other”  compartment. 


Appendix  III.  PBPK  Plot  Simulations 


Plot  V:  High  dosed  male  medaka  using  high  dose  male  concentration  ratios  as  partition  coefficients,  Vmaxc  set  at  0.128 


Appendix  III.  PBPK  Plot  Simulations 


VI:  Low  dosed  male  medaka  using  high  dose  male  concentraton  ratios  as  partition 
coefficients,  Vmaxc  set  at  0.128  (unacceptable,  CL  and  CG  far  outside  lsd  range). 
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/II:  Simulation  of  low  dosed  male  medaka,  using  low  dose  male  concentration  ratios  as  partitions,  and  Vmaxc 
0.128  (an  acceptable  result). 


Appendix  III.  PBPK  Plot  Simulations 


Plot  VIII:  Simulation  of  high  dosed  male  medaka,  using  low  dose  male  concentration  ratios  as  partitions,  and  Vmaxc  set  at 
0.128.  (an  unacceptable  result  due  to  the  fact  that  CV  simulated  concentrations  were  far  outside  the  lsd  range). 


Appendix  III.  PBPK  Plot  Simulations 


Plot  IX:  Simulation  of  high  dose  male  medaka,  using  low  dosed  male  concentration  ratios  as  partitions, 
and  Vmaxc  set  at  0.0128  (unacceptable  due  to  CV  and  CM  both  falling  outside  their  lsd  ranges) 


Appendix  II 


0.0128  (acceptable). 


Appendix  IV.  Statistical  values  for  histopathology  endpoints  of  Test  108-002  (highlighted  values  are  significant  at  p  =  0.05) 
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Appendix  IV.  Statistical  values  for  histopathology  endpoints  of  Test  108-002  (highlighted  values  are  significant  at  p  =  0.05) 


IV-3 


Appendix  IV.  Statistical  values  for  histopathology  endpoints  of  Test  108-002  (highlighted  values  are  significant  at  p  =  0.05) 
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